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Abstract 



The observation of the exotic quarkonium state Zc(3900) by the BESIII and Belle collaborations supports the concept of hadronic 
molecules. Charmonium states interpreted as such molecules would be bound states of heavy particles with small binding energies, 
f^ This motivates their description using an effective theory with contact interactions. In particular, we focus on the electromagnetic 
T-H structure of the charged state Zc(3900). Using first experimental results concerning spin and parity, we interpret it as an 5 -wave 
^^ molecule and calculate the form factors as well as charge and magnetic radii up to next-to-leading order. We also present first 
^^ numerical estimations of some of these observables at leading order 






Introduction 



In 2013 the BESIII and Belle collaborations reported the ob- 
I— I servation of a charged state in the charmonium sector called 
^ Zc(3900) 1 1 , 2 1 (this state was also seen at CLEG Lll). The 
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measured mass and width are 



Mz = (3899.0 ± 8.5) MeV Tz = (46 ± 30) MeV . (1) 



Soon after its discovery it was proposed that these states are S - 
Cn wave hadronic molecules |l4l|5l- Later a further analysis in this 
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framework and different approaches for the substructure like 
tetraquarks were presented in Ref. |6|. In this work, we will 
assume that the Zc is indeed a hadronic molecule. A comparison 
with combinations of charm mesons shows that this states lies 
close to a specific threshold and thus favors this interpretation. 
Denoting the wave function with the respective particle name. 



C^ one can identify their main component to be 141161 



X 



Zc = 



V2 



{pD* + DD") . 



(2) 



First, the constituents are charm mesons whose masses are far 
above 1 GeV and second, their binding energy is rather small, 
i.e. below the pion mass m„ or at least of that order Hence, 
for the description of these states one can use in a first approx- 
imation a pionless effective field theory where the constituents 
are treated as non-relativistic, point-like particles which only 
interact via contact interactions. Such an short-range effective 
theory is called EFIiij:) and is a commonly used tool in nuclear 
physics |7-10|. Its expansion parameter is Q/m„, where Q is 
the scale of the modulus p of the internal momenta p of the 
involved particles which are - according to the small binding 
energy of the molecule, \B\ < 30 MeV - below the pion mass 
m„. We already remark here that the characteristics of the Zc 
are such that EFT(:^) is barely applicable, but it still allows for 



some first insights into the electromagnetic structure of this ex- 
otic particle. This can be made more quantitative by explicitly 
including pions in the underlying EFT. 

For the charm molecule Z^ the analysis of its decay chan- 
nel shows that the quantum numbers /(7^) = l(l^) are favored 
HI El, which agrees with the 5 -wave hypothesis. This fact al- 
lows one to relate the structure of the Zc state to that of the 
deuteron treated in EFT(:}<-) as a baryonic molecule of nucleons 
in a relative S -wave with total spin 7=1 fTr-T?!. Besides this 
analogy, there is also a conceptual connection to the so-called 
Halo EFT developed in Refs. lfT6l [TTll and in particular to its 
application in Ref. lTSl . 



2. Formalism 

Following Refs. ifTTl [TSl . one can write the effective La- 
grangian up to next-to-leading order (NLG) for the positively 
charged Zc as 
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Figure 1: Full Z propagator (thick solid line) as sum over charm meson loops. 
The double solid line represents the bare Z propagator, the single solid line 
stands for D* mesons, the dashed one for D mesons and the dots represent 
higher order terms. 



Here, Z is an auxiliary field describing the molecular state 
Zc(3900), g is the coupling constant for ZDD* interactions, A 
is a constant and rf - ±1 is a phase whose sign depends on the 
effective range. Minimal substitution with electric charge e in- 
duces an Aq interaction with the constituents and the molecule 
itself, leading to the Feynman diagrams which contribute to the 
charge form factor. The magnetic form factor instead comes 
from the terms in the last line of Eq. ([3]), where B = V x A 
is the magnetic field and U is the S0(3) generalization of the 
PauU vector, with t/,- being the generators of the rotation group. 
Those terms are proportional to the magnetic moment fi of the 
respective charged particle and describe purely magnetic inter- 
actions of the spatial part of the vector potential A,- with both 
the charged D mesons and Z^. Other contributions to the La- 
grangian, for instance those relevant for the quadrupole form 
factor or effects due to the anomalous magnetic moment of neu- 
tral D mesons, are at least one order higher and thus not shown. 
Also the explicit spin structure of the Lagrangian can be omit- 
ted since the D meson is a scalar, i.e. the Z and D* are always 
in the same spin state. 

The auxiliary field Z is not dynamical but due to its coupling 
to the constituents it is dressed by D meson loops, see Fig.[T] 
This leads to a full Z propagator 



iSz - 



{sr-^ 



(4) 



in terms of the bare propagator iS^^ and the self-energy E. To 
avoid an unnatural scaling of the scattering length (similar to 
that discussed in Ref. IIT8II ). one calculates the self-energy in 
the PDS scheme introduced in Refs. [,19n20J. In terms of the 
PDS scale Apos and the mean reduced mass toz = j( 



Mo+M,4- 
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2 *- M++M,o 



) of the molecule, one obtains 



S = - 



(^zg 



i ^j2cjz \pq- —\ + is + ApDs 



(5) 



From the full propagator of the Z field at NLO, 



iSz(po,p) 



(^zg^ 



(^zg^ 



+ ApDs + iJ2uJz\po 
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cozg^ y"'' 2Mz 
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+ IS 



,(6) 



one can deduce the scattering length a and the effective 
range ro of the Zc by matching the scattering amplitude 



-iT = (-ig)^iSz(E = k^/{2a)z),0) to the effective range ex- 
pansion (ERE) 
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At NLO the following relations are obtained: 
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(7) 

(8) 
(9) 



In addition this result fixes the phase introduced in the La- 
grangian Eq. ([3]l. To get a positive value for the effective range, 
the phase j] must be equal to minus one. 

Furthermore, we calculate the wave function renormalization 
constant Z defined as the residue of the pole of the propagator 
in Eq. (|6]l. At NLO one finds 

1 
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where the Zc binding momentum y - sgn(B) y]2(jL>z\B\ was 
introduced whose definition is chosen such that one takes care 
of both, bound and virtual states. 

The non-relativistic matrix elements of the electromagnetic 
current y^ can be expanded in terms of form factors. For a 
vector particle state |p, e,) with three momentum p and linear 
polarization vector e,, one can write in the Breit-frame lfTTl [151 



<p',ej|7"|p,e,) = -ie[Gc{q)6ij 
1 
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(11) 
(12) 



Here, q = p' - p is the three -momentum transfer with q^q'^ = 
q^ - cf- and n is the number of space-time dimensions (here 
n - 4). The charge {Gc), magnetic (G,„) and quadrupole (Gq) 
form factors are normalized in the following way 121]: 



Gc{q 
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(13) 



in terms of the magnetic moment [i of the particle and its 
quadrupole moment k. 



According to Eq. (Ill, the charge form factor is determined 
by the diagrams shown in Fig.l2]in which the incoming photons 
are Aq photons and where the initial and final spin of the Zc are 
equal. The charge form factor at NLO of the Zc is then given by 



G^^Kq) 



1 



1 -yo 

M.+y 



qujz 

i 



arctan 



arctan 



qb)z 



qoz 



\2M,+y 



qi^z 

2M+y 



- yo 



(14) 




Figure 2: Diagrams contributing to tlie cliarge and magnetic form factor of the 
Zr at NLO. While the left diagrams in parenthesis are LO, the right diagram 
is the NLO contribution. The thick solid line represents the full Z propagator, 
the thin solid line stands for D* and the dashed line for D mesons. The incom- 
ing photons are indicated by wavy lines and the application of wave function 
renormalization factors is indicated by the crosses. 



3. Results 

Except for the masses, widths and liJ'') there are no exper- 
imental data available for the Zc(3900), hence one can predict 
observables only at LO. Furthermore, the leading order mag- 
netic form factor and the magnetic radius are proportional to 
the magnetic moments of the constituent mesons. As these are 
experimentally unknown both are not accessible even at LO. 
Thus only the charge form factor and the charge radius at LO 
can be predicted. 

Using the charm meson masses from Ref. Il22l the following 
(averaged) binding energy B can be obtained 
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This is indeed correctly normalized to one. A second observ- 
able besides the form factor itself is the expectation value of the 
squared charge radius {rl). It can be found from the low-energy 
expansion of the form factor, Gdq) = 1 - {rl)q^ 16 + 0{q*). At 
NLO one finds 



<r?) 
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(15) 



If the incoming photons in Fig.l2]are A, photons, one can choose 
the polarizations in Eq. (fT2]i and calculate the amplitude of the 
corresponding diagrams to get the magnetic form factor G,„ of 
the Ze. At NLO it is given by 



■^(NLO) 
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(16) 

So the difference to the charge form factor are the magnetic 
moments in front of each term. With the normalization con- 
dition in Eq. ( [T3] l. it is possible to identify a relation between 
the magnetic moment of the molecule Z^ and that of its charged 
constituents D^ and D*^ 



f^z = 7(y"+ +ju*+) 



(17) 



In addition, one can determine the expectation value of the 
squared magnetic radius (r^,) by expanding the magnetic form 
factor up to second order in q. At NLO one finds 
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(18) 



According to Eq. (Ill, the quadrupole form factor is determined 
by A() photon interactions in Fig. l2]with a spin change between 
the initial and the final state. The Aq interaction itself cannot 
induce such a spin change because Aq is just the time-like com- 
ponent of the 4- vector potential A^. Thus one needs additional 
operators in the Lagrangian which project on specific polariza- 
tion states e,. According to Ref. 1 1 1 1 such operators first appear 
at NLO and thus the quadrupole form factor Gq{q) vanishes at 
leading order 



B = - [(M+ + M,o - Mz) + (Mo + M,+ - Mz)] 
= (-23.13 + 6.01) MeV. 



(19) 



The negative value contradicts the charm meson bound state so- 
lution and favors a resonance. Neglecting the imaginary part of 
the binding momentum y, we treat the Zc(3900) as a a virtual 
state with negative scattering length. In leading order, the elec- 
tromagnetic properties are independent of the sign of y. Our 
EFT expansion parameter is Q/m„ and since the inverse scat- 
tering length 1 /a should be a low-energy scale proportional to 
Q to provide a proper counting scheme, one concludes that |a| 
must be larger than 1 /m„ a; 1 .4 fm. This is not the case for the 
scattering length of the Z^ given below. However, one can use 
the predictions for the Z^. as first, rough model estimations. 

The results for the binding momentum, the scattering length 
and charge radius of the Z^ are 



y = (-211.6 + 27.5) MeV ^ 
a = (-0.93 ± 0.12) fm, 
(r^) = (0.11 ±0.03)fm^ 



(20) 



where the uncertainties are entirely from the one in the binding 



energy, Eq. ( 19 1. Due to the large expansion parameter, we re- 
frain from estimating the errors due to higher orders. Also, note 
that the charge radius in Eq. pO] ) should be added in quadrature 
the charge radii of the constituent D^*' mesons which are treated 
as pointlike in our theory. 

The charge form factor of the Z^ at LO is shown in Fig. [3| 
Furthermore, there is a the correlation between the leading or- 
der charge radius and the binding energy as shown in Fig. |4] 

4. Summary 

In this work the electromagnetic structure of the exotic me- 
son Zc(3900) interpreted as 5 -wave charm meson molecule 
with spin J - 1 was investigated in the framework of EFT(.5(-). 
Based on the description of the deuteron in this eff'ective theory, 
the charge, magnetic and quadrupole form factor were formally 
analyzed up to NLO. In particular, the charge and magnetic ra- 
dius were obtained. Due to the poor experimental data on the 
Zc state only the LO charge radius and charge form factor could 
be estimated. Furthermore, we derived a relation between the 
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Figure 3: LO charge form factor gJ. (q) of the Z£.(3900). The error band 
includes only the uncertainty in the binding energy and is depicted in green. 

magnetic moment of the molecule and of its constituents, see 
Eq. ^Vl\ . Given the experimental mass of the Zc(3900), the 
binding momentum of the 2^(3900) is a bit too large to allow for 
a quantitative description of its properties in EFT(:5<-). Thus our 
results should be viewed as as first, rough model estimations 
of the electromagnetic properties of the Zc(3900). To sharpen 
these predictions, one has to go beyond the pionless EFT and 
include explicit pions. 
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Figure 4: Correlation between the expectation value of the squared charge ra- 
dius (r~) and the modulus of the binding energy B of the Zc(3900) at LO. The 
congelation is shown for binding energies within the eiTors of Eq. jl9^ . The 
central value Eq. l |20) is represented by the cross and the error band is depicted 
in green. 
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